Population phylogenomic analysis of mitochondrial DNA in wild boars and domestic pigs revealed multiple domestication events in East Asia

Background: Previously reported evidence indicates that pigs were independently domesticated in multiple places throughout the world. However, a detailed picture of the origin and dispersal of domestic pigs in East Asia has not yet been reported. Results:

The evolution of YidC/Oxa/Alb3 family in the three domains of life: a phylogenomic analysis

Background: YidC/Oxa/Alb3 family includes a group of conserved translocases that are essential for protein insertion into inner membranes of bacteria and mitochondria, and thylakoid membranes of chloroplasts. Because mitochondria and chloroplasts are of b

Global phylogenomic analysis of nonencapsulated Streptococcus pneumoniae reveals a deep-branching classic lineage that is distinct from multiple sporadic lineages.

The surrounding capsule of Streptococcus pneumoniae has been identified as a major virulence factor and is targeted by pneumococcal conjugate vaccines (PCV). However, nonencapsulated Streptococcus pneumoniae (Non-Ec-Sp) have also been isolated globally, mainly in carriage studies. It is unknown if Non-Ec-Sp evolve sporadically, if they have high antibiotic non-susceptiblity rates and a unique, specific gene content. Here, whole genome sequencing of 131 Non-Ec-Sp isolates sourced from 17 different locations around the world was performed. Results revealed a deep-branching classic lineage that is distinct from multiple sporadic lineages. The sporadic lineages clustered with a previously sequenced, global collection of encapsulated S. pneumoniae (Ec-Sp) isolates while the classic lineage is comprised mainly of the frequently identified multi-locus sequences types ST344 (n=39) and ST448 (n=40). All ST344 and nine ST448 isolates had high non-susceptiblity rates to β-lactams and other antimicrobials. Analysis of the accessory genome reveals that the classic Non-Ec-Sp contained an increased number of mobile elements, than Ec-Sp and sporadic Non-Ec-Sp. Performing adherence assays to human epithelial cells for selected classic and sporadic Non-Ec-Sp revealed that the presence of a integrative conjugative element (ICE) results in increased adherence to human epithelial cells (P=0.005). In contrast, sporadic Non-Ec-Sp lacking the ICE had greater growth in vitro possibly resulting in improved fitness. In conclusion, Non-Ec-Sp isolates from the classic lineage have evolved separately. They have spread globally, are well adapted to nasopharyngeal carriage and are able to coexist with Ec-Sp. Due to continued use of pneumococcal conjugate vaccines, Non-Ec-Sp may become more prevalent.

Phylogenomic analysis of natural products biosynthetic gene clusters allows discovery of arseno-organic metabolites in model streptomycetes

We are indebted with Marnix Medema, Paul Straight and Sean Rovito, for useful discussions and critical reading of the manuscript, as well as with Alicia Chagolla and Yolanda Rodriguez of the MS Service of Unidad Irapuato, Cinvestav, and Araceli Fernandez for technical support in high-performance computing. This work was funded by Conacyt Mexico (grants No. 179290 and 177568) and FINNOVA Mexico (grant No. 214716) to FBG. PCM was funded by Conacyt scholarship (No. 28830) and a Cinvestav posdoctoral fellowship. JF and JFK acknowledge funding from the College of Physical Sciences, University of Aberdeen, UK.

Epidemiology and comparative analysis of Yersinia in Ireland

Yersiniosis is an acute or chronic enteric zoonosis caused by enteropathogenic Yersinia species. Although yersiniosis is predominantly associated with gastroenteric forms of infection, extraintestinal forms are often reported from the elderly or patients with predisposing factors. Yersiniosis is often reported in countries with cold and mild climates (Northern and Central Europe, New Zealand and North of Russian Federation). The Irish Health Protection Surveillance Centre (HPSC) currently records only 3-7 notified cases of yersiniosis per year. At the same time pathogenic Yersinia enterocolitica is recovered from pigs (main source of pathogenic Y. enterocolitica) at the levels similar to that observed in Yersinia endemic countries. Introduction of Yersinia selective culture procedures may increase Yersinia isolation rates. To establish whether the small number of notifications of human disease was an underestimate due to lack of specific selective culture for Yersinia we carried out a prospective culture study of faecal samples from outpatients with diarrhoea, with additional culture of appendix and throat swabs. Higher levels of anti-Yersinia seroprevalence than yersiniosis notification rates in endemic countries suggests that most yersiniosis cases are unrecognised by culture. Subsequently, in addition to a prospective culture study of clinical specimens, we carried out serological screening of Irish blood donors and environmental screening of human sewage. Pathogenic Yersinia strains were not isolated from 1,189 faeces samples, nor from 297 throat swabs, or 23 appendix swabs. This suggested that current low notification rates in Ireland are not due to the lack of specific Yersinia culture procedures. Molecular screening detected a wider variety of Y. enterocolitica-specific targets in pig slurry than in human sewage. A serological survey for antibodies against Yersinia YOP (Yersinia Outer Proteins) proteins in Irish blood donors found antibodies in 25%, with an age-related trend to increased seropositivity, compatible with the hypothesis that yersiniosis may have been more prevalent in Ireland in the recent past. Y. enterocolitica is a heterogeneous group of microorganisms that comprises strains with different degree of pathogenicity. Although non-pathogenic Y. enterocolitica lack conventional virulence factors, these strains can be isolated from patients with diarrhoea. Insecticidal Toxin Complex (ITC) and Cytolethal Distending Toxins can potentially contribute to the virulence of non-pathogenic Y. enterocolitica in the absence of other virulence factors. We compared distribution of ITC and CDT loci among pathogenic and non-pathogenic Y. enterocolitica. Additionally, to demonstrate potential pathogenicity of non-pathogenic Y. enterocolitica we compared their virulence towards Galleria mellonella larvae (a non-mammalian model of human bacterial infections) with the virulence of highly and mildly pathogenic Y. enterocolitica strains. Surprisingly, virulence of pathogenic and non-pathogenic Y. enterocolitica in Galleria mellonella larvae observed at 37°C did not correlate with their pathogenic potential towards humans. Comparative phylogenomic analysis detects predicted coding sequences (CDSs) that define host-pathogen interactions and hence providing insights into molecular evolution of bacterial virulence. Comparative phylogenomic analysis of microarray data generated in Y. enterocolitica strains isolated in the Great Britain from humans with diarrhoea and domestic animals revealed high genetic heterogeneity of these species. Because of the extensive human, animal and food exchanges between the UK and Ireland the objective of this study was to gain further insight into genetic heterogeneity and relationships among clinical and non-clinical Y. enterocolitica strains of various pathogenic potential isolated in Ireland and Great Britain. No evidence of direct transfer of strains between the two countries was found.

Missed, not missing: phylogenomic evidence for the existence of Avian FoxP3

The Forkhead box transcription factor FoxP3 is pivotal to the development and function of regulatory T cells (Tregs), which make a major contribution to peripheral tolerance. FoxP3 is believed to perform a regulatory role in all the vertebrate species in which it has been detected. The prevailing view is that FoxP3 is absent in birds and that avian Tregs rely on alternative developmental and suppressive pathways. Prompted by the automated annotation of foxp3 in the ground tit (Parus humilis) genome, we have questioned this assumption. Our analysis of all available avian genomes has revealed that the foxp3 locus is missing, incomplete or of poor quality in the relevant genomic assemblies for nearly all avian species. Nevertheless, in two species, the peregrine falcon (Falco peregrinus) and the saker falcon (F. cherrug), there is compelling evidence for the existence of exons showing synteny with foxp3 in the ground tit. A broader phylogenomic analysis has shown that FoxP3 sequences from these three species are similar to crocodilian sequences, the closest living relatives of birds. In both birds and crocodilians, we have also identified a highly proline-enriched region at the N terminus of FoxP3, a region previously identified only in mammals.

Insect mitochondrial genomics : implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

蔷薇类的分子系统学研究—兼论系统发育基因组学的原理和方法(以SET基因超家族为例)

分子系统发育分析的主要任务包括：（1）帮助建立生命之树（tree of life）;（2）追踪基因和基因家族（gene family）的起源和进化, 以获知基因在进化过程中的功能分化和伴随发生的重要分子事件（key molecular events）和形态性状的关键创新（key innovation）。这两个方面在本研究中都有所涉及。对于前者，选用植物线粒体matR基因重建被子植物蔷薇类群的系统发育关系;对于后者，则以SET基因超家族为例，探讨其在真核生物中的进化分类以及与功能多样性的关系。 I 蔷薇类的分子系统学 蔷薇类（rosids）是基于分子数据建立的被子植物的主要分支之一，包含13个目，大约三分之一的被子植物物种。两个主要蔷薇类内部分支是豆类fabids（包含7个目）和锦葵类malvids（包含3个目）。然而，这两个分支内部，以及这两个分支与蔷薇类基部类群，包括牻牛儿苗目（Geraniales）、桃金娘目（Myrtales）和流苏子目（Crossosomatales）之间的关系大多是不清楚的。本研究中，我们选取174个物种来代表72个蔷薇类（rosids）的科，利用两个数据集，即线粒体matR单基因数据集和包括线粒体matR基因、两个质体基因(rbcL、 atpB)和一个核基因(18S rDNA) 的4基因数据集，重建蔷薇类在科以上分类阶元水平的系统发育关系。同时，还对线粒体matR基因的进化特征和用于大尺度系统发育分析的适合度和潜力进行了评价。 线粒体matR单基因数据支持malvids和大多数蔷薇类目的单系性质，然而，豆类（fabids）成员没有形成一个分支，其COM亚支，包括卫矛目（Celastrales）、酢浆草目（Oxalidales）、金虎尾目（Malpighiales）和蒜树科（Huaceae），分辨为锦葵类（malvids）的姐妹群。这个关系在最近根据花结构特征曾被提出过，但从未在之前的分子系统发育分析中得到分辨。4基因数据集支持首先是牻牛儿苗目（Geraniales），接着是桃金娘目（Myrtales）作为蔷薇类（rosids）的最基部的分支;流苏子目（Crossosomatales）是锦葵类（malvids）姐妹群，以及蔷薇类（rosids）的核心部分包括豆类（fabids），锦葵类（malvids）和流苏子目（Crossosomatales）。线粒体matR基因的进化特征分析显示，与两个叶绿体基因(rbcL 和atpB)比较，同义替代速率约是它们的1/4，而非同义替代速率接近于自身的同义替代速率，表明matR 基因具有松弛的选择压力。线粒体matR基因相对慢速的进化使非同源相似（homoplasious）突变减少，提高了系统发育信息的质量，同时，松弛的选择压力使非同义替代数量增加，弥补了慢速进化导致的系统发育信息数量不足的缺陷，这两个方面的结合使线粒体matR基因非常适用于被子植物在科以上水平的系统发育研究。 II SET基因超家族的系统发育基因组学分析 SET基因超家族基因编码含有SET结构域的蛋白，在真核生物中，SET-domain蛋白一般是多结构域（multi-domain）的。SET-domain蛋白具有对组蛋白H3和H4的N末端尾部进行赖氨酸残基甲基化修饰的酶活性;从异染色质形成到基因转录，甲基化的组蛋白广泛影响染色质水平的基因调控。依据SET结构域一级序列的相似性和结构域组织（domain architecture）特征，目前，SET-domain基因超家族被划分为4-7个家族。由于这些划分或者使用动物或者使用植物SET基因，只有少数其它类群的物种加入分析，因此这样的划分可能是不完整的。本研究采用系统发育基 因组学方法（phylogenomic approach），在真核生物范围内广泛取样，期望获得相对完整的SET-domain基因家族的 进化分类方案，在此基础上加深理解SET-domain基因的进化机制和功能多样性。 在提取了17个物种，代表5个真核超群的SET蛋白序列基础上，系统发育分析结合“结构域组织特征”鉴别了9个SET基因家族，其中一个是新的SET基因家族。以前的SET8和Class VI家族,及SMYD和SUV4-20家族分别合并为一个家族。大部分家族在进化过程中发生了2次以上的基因重复事件，通过获得不同的结构域产生具有不同功能的新基因。一个SET基因家族在进化过程中推测发生了从脊椎动物祖先向盘基网柄菌（Dictyostelium discoideum）的水平基因转移。

Comparative Genomics of Early-Diverging Brucella Strains Reveals a Novel Lipopolysaccharide Biosynthesis Pathway

Brucella species are Gram-negative bacteria that infect mammals. Recently, two unusual strains (Brucella inopinata BO1(T) and B. inopinata-like BO2) have been isolated from human patients, and their similarity to some atypical brucellae isolated from Australian native rodent species was noted. Here we present a phylogenomic analysis of the draft genome sequences of BO1(T) and BO2 and of the Australian rodent strains 83-13 and NF2653 that shows that they form two groups well separated from the other sequenced Brucella spp. Several important differences were noted. Both BO1(T) and BO2 did not agglutinate significantly when live or inactivated cells were exposed to monospecific A and Mantisera against O-side chain sugars composed of N-formyl-perosamine. While BO1(T) maintained the genes required to synthesize a typical Brucella O-antigen, BO2 lacked many of these genes but still produced a smooth LPS (lipopolysaccharide). Most missing genes were found in the wbk region involved in O-antigen synthesis in classic smooth Brucella spp. In their place, BO2 carries four genes that other bacteria use for making a rhamnose-based O-antigen. Electrophoretic, immunoblot, and chemical analyses showed that BO2 carries an antigenically different O-antigen made of repeating hexose-rich oligosaccharide units that made the LPS water-soluble, which contrasts with the homopolymeric O-antigen of other smooth brucellae that have a phenol-soluble LPS. The results demonstrate the existence of a group of early-diverging brucellae with traits that depart significantly from those of the Brucella species described thus far. IMPORTANCE This report examines differences between genomes from four new Brucella strains and those from the classic Brucella spp. Our results show that the four new strains are outliers with respect to the previously known Brucella strains and yet are part of the genus, forming two new clades. The analysis revealed important information about the evolution and survival mechanisms of Brucella species, helping reshape our knowledge of this important zoonotic pathogen. One discovery of special importance is that one of the strains, BO2, produces an O-antigen distinct from any that has been seen in any other Brucella isolates to date.

Dissimilarities in alignment-free methods for phylogenetic analysis based on genomes

Whole genome sequences are generally accepted as excellent tools for studying evolutionary relationships. Due to the problems caused by the uncertainty in alignment, existing tools for phylogenetic analysis based on multiple alignments could not be directly applied to the whole-genome comparison and phylogenomic studies. There has been a growing interest in alignment-free methods for phylogenetic analysis using complete genome data. The “distances” used in these alignment-free methods are not proper distance metrics in the strict mathematical sense. In this study, we first review them in a more general frame — dissimilarity. Then we propose some new dissimilarities for phylogenetic analysis. Last three genome datasets are employed to evaluate these dissimilarities from a biological point of view.

Phylogenomic study of the subfamily Caprinae by cross-species chromosome painting with Chinese muntjac paints

Chromosomal homologies have been established between the Chinese muntjac (Muntiacus reevesi, MRE, 2n = 46) and five ovine species: wild goat (Capra aegagrus, CAE, 2n = 60), argall (Ovis ammon, OAM, 2n = 56), snow sheep (Ovis nivicola, ONI, 2n = 52), red goral (Naemorhedus cranbrooki, NCR, 2n = 56) and Sumatra serow (Capricornis sumatraensis, CSU, 2n = 48) by chromosome painting with a set of chromosome-specific probes of the Chinese muntjac. In total, twenty-two Chinese muntjac autosomal painting probes detected thirty-five homologous segments in the genome of each species. The chromosome X probe hybridized to the whole X chromosomes of all ovine species while the chromosome Y probe gave no signal. Our results demonstrate that almost all homologous segments defined by comparative painting show a high degree of conservation in G-banding patterns and that each speciation event is accompanied by specific chromosomal rearrangements. The combined analysis of our results and previous cytogenetic and molecular systematic results enables us to map the chromosomal rearrangements onto a phylogenetic tree, thus providing new insights into the karyotypic evolution of these species.

Raman microscopy of formamide-intercalated kaolinites treated by controlled-rate thermal analysis technology

Raman spectroscopy of formamide-intercalated kaolinites treated using controlled-rate thermal analysis technology (CRTA), allowing the separation of adsorbed formamide from intercalated formamide in formamide-intercalated kaolinites, is reported. The Raman spectra of the CRTA-treated formamide-intercalated kaolinites are significantly different from those of the intercalated kaolinites, which display a combination of both intercalated and adsorbed formamide. An intense band is observed at 3629 cm-1, attributed to the inner surface hydroxyls hydrogen bonded to the formamide. Broad bands are observed at 3600 and 3639 cm-1, assigned to the inner surface hydroxyls, which are hydrogen bonded to the adsorbed water molecules. The hydroxyl-stretching band of the inner hydroxyl is observed at 3621 cm-1 in the Raman spectra of the CRTA-treated formamide-intercalated kaolinites. The results of thermal analysis show that the amount of intercalated formamide between the kaolinite layers is independent of the presence of water. Significant differences are observed in the CO stretching region between the adsorbed and intercalated formamide.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.